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The transport properties of carriers in semiconducting graphene nanoribbons are studied by com- 
paring the effects of phonon, impurity, and line-edge roughness scattering. It is found that scattering 
from impurities located at the surface of nanoribbons, and from acoustic phonons are as important 
as line edge roughness scattering. The relative importance of these scattering mechanisms varies 
with the temperature, Fermi level location, and the width of the ribbons. Based on the analysis, 
strategies for improvement of low- field mobility are described. 
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I. INTRODUCTION 



Since the isolation of 2-dimensional (2D) graphene 
sheets in 2004 [1, 2, 3], there has been substantial interest 
in patterning them into quasi- ID graphene nanoribbons 
(GNRs) [4]. This interest is motivated by the ability to 
tune the bandgap of GNRs by lithographic processes, and 
the consequent freedom it is expected to afford in litho- 
graphic bandgap-engineered electronic devices. GNRs as 
thin as ~ 3 nm have been achieved, exhibiting bandgaps 
approaching 0.3 - 0.4 eV, which makes them attractive for 
many low-power device applications [5]. The measured 
carrier mobilities in the first reported ultrathin GNRs are 
reported to be much lower [6] than in corresponding 2D 
graphene sheets [7]. For improvements in the transport 
properties, it is essential to identify the major sources of 
carrier scattering in thin GNRs. To that end, this work 
presents a comprehensive analysis of various scattering 
mechanisms in GNRs. Analytical scattering rates for 
phonon, impurity, and line-edge roughness (LER) scat- 
tering are derived taking the GNR wavefunctions into 
account. It is shown that carrier mobilities in GNRs are 
not necessarily limited by edge roughness scattering; sim- 
ilar to the case in 2D graphene, surface impurity [8] and 
phonon scattering play important roles, depending on the 
carrier concentrations, widths, and temperature. The re- 
sults of this work are expected to provide useful strategies 
towards improvements in carrier transport properties in 
GNRs for various applications. 

The paper is structured in the following fashion. Sec- 
tion II describes the theoretical formalism used for the 
scattering rate and mobility calculations, with a discus- 
sion of the electronic structure, scattering matrix ele- 
ments, and screening in GNRs. In section III, individual 
scattering rates due to acoustic and optical phonons, line 
edge roughness (LER), and bulk and surface impurities 
are derived. Finally, in section IV, the resulting mobility 
in GNRs is calculated, with a discussion of the relative 



importance of the various scattering mechanisms. The 
effect of temperature, location of Fermi level, and GNR 
width on the transport properties is evaluated in that 
section. 



II. THEORETICAL FORMALISM 

Under the application of a small electric field along 
the GNR axis, the carrier distribution function in the 
relaxation time approximation is given by 
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where r _1 is the scattering rate in the diffusive limit, v g 
is the group velocity of carriers, T is the electric field, e 
is the electron charge, and /o is the equilibrium Fermi- 
Dirac distribution. If the scattering rate is known, the 
current flowing through the GNR may be evaluated by 
the relation 
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where e is the electron charge, L is the GNR length, 
g s = 2 is the spin degeneracy and g v — 1 is the valley 
degeneracy of semiconducting GNRs. Noting that the 
electric field is T — V/L where V is the applied voltage, 
the ID conductivity (in units of S-m) is given by 
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where h = 2irh is Planck's constant. The mobility is 
calculated as /i = (Jid/^^id^ where nuj is the one- 
dimensional carrier density. This formalism allows for 
the evaluation of scattering rates, mobility, and conduc- 
tivity for a general Fermi level location, which can be 
assumed to be tuned capacitively through a gate voltage 
in experiments, at any temperature. 



A. Electronic Structure 

Semiconducting armchair GNRs with lengths L(» 
W) along the y (x) directions are considered in this 
work. The hard-wall boundary conditions at the edges 



lead to a bandstructure £(k n ,k y ) = hvpJk^ + k^ where 

vp ~10 8 cm/s is the Fermi velocity. The allowed trans- 
verse wavevector for semiconducting armchair GNRs are 
quantized to values k n = ±n7r/3W, depending on the 
width W; here n = ±1, ±2, ±4, ±5, ±7, ±8, .... The re- 
sulting bandgap is £ g = 2ttHvf/SW : which is ~ 1.38 /W 
eV with W expressed in nm. We note that the group 
velocity of carriers in a GNR is linked to the density of 
states (DOS) pgnr by v g (£ ) = 2/7rhpGNR(£), where the 
DOS is p GNR (S) = (2/7rhv F ) x (£/^£ 2 - (£ g /2) 2 ) for 
the first subband. The effective ID carrier density in the 
GNR is given by n 1D = J °° pGNR(£)fo(£)d£. 



B. Matrix Elements and Scattering Rates 

The electronic wavefunctions of carriers in semicon- 
ducting GNRs may be derived from an admixture of 
states in the JC and K! valleys of the underlying 2D 
graphene bandstructure [9]. For armchair semiconduct- 
ing GNRs, the valley degeneracy is removed. Assuming 
a length L and a width W ", the electronic wavefunctions 
in the nth subband in the Ket-notation are written as 
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where the tilde sign indicates the state in the K! val- 
ley. The resulting energy dispersion (the £ — k relation) 

of the nth subband is £(k nj k y ) = hvpJk^ + fc 2 where 

vf ~ 10 8 cm/s is the Fermi velocity in graphene. The 
projections of these states on to the real space gives the 
wavefunctions 
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where r = (x, y) is the vector in the x-y plane. AK = 
47r/3a is the distance between two Dirac points in the 
fc-space for the underlying 2D graphene bandstructure, 
a being the lattice constant of graphene. N^- — k n is 
the electron wave vector in the confined x-direction for 
the nth subband states. n is the angle between k y and 
K m @n — tan _1 (/c 2y //c n ). Since the carrier states in the 
/c-space are close to the Dirac point, k n « AK. 

For a perturbation of the form V(x, y, z), the pertuba- 



tion matrix element is given by 

(<p'\V(x,y,z)\tp) = \{(k' y ,k n ,\V\k n ,k v ) + 
\ky> k n i | V | k n , ky ) > — 

\{(k' y X'\V\k n ,ky) + (k' y ,k n/ \V\kZ,ky)}. (6) 

For finding the mobility and conductance, the scatter- 
ing rates due to various scattering mechanisms are eval- 
uated using Fermi's golden rule, in the form 

2-7T 

5(k,k') = T K^in^y^)k)l 2 ^k-4±M, (7) 

where the delta function ensures energy conservation for 
both elastic (uo = 0) inelastic (uo ^ 0) scattering pro- 
cesses. The ensemble scattering rate that contributes to 
the conductivity and mobility is then evaluated for each 
scattering mechanism as 



- = ^S(k,k / )(l-cosa), 



(8) 



a is the angle between k and k 7 , and the summation runs 
over all available final states. 



C. Screening 

Screening of the scattering potential by free-carriers 
in the GNRs modify the scattering rates. To take this 
many-body effect into account, the following procedure 
has been used (see [10]). The static screening is cal- 
culated using the random phase approximation (RPA). 
In the k— space, the relation between screened potential 
and unscreened potential is determined by V uns (q n ) = 
E m e (<7m, qn)V scr (q m )' The screening matrix is given by 

e 2 

e(qm,q n ,q y ) = ^ m gn + o 3 r (qm,qn,q y )£(qn,q y ), (9) 
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where eo is the permittivity of vacuum, k is the average 
of the dielectric constants of the regions between which 
the GNR is sandwiched, and S qrn q n is the Kronecker sign. 
!FQ} is a measure of the contribution to screening from 
coupled transverse modes, given by 

3 r (qm,qn,q y ) = ^2 / / K (\q y (x - x')\) x 

cos (q m x) cos (q n x f )dxdx f \ (10) 

where Kq(...) is the zero-order modified Bessel function. 
C(q n ,q y ) is given by 
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which is a sum over Lindhard functions C nn > between 
two subbands that have a wave vector difference q n in 



the transverse (x— direction). This sum includes both in- 
tersubband and interband (conduction and valence band) 
contributions. The Lindhard function is given by 
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where s = 1 is for conduction subbands and s' = — 1 is 
for valence subbands. g s = 2 is the spin degeneracy, and 
/n(--) is the Fermi-Dirac distribution function in the nth 
subband. If \q m — q n \ = Itt/W and I is an odd number, 
^(Qm^Qn^Qy) = 0. This makes the screening calculation 
very simple for the lowest subbands. The screening of 
the first subband is e scr = e(0,0,g y ), which is given by 
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which is equal to 1 if screening is neglected. In the 
long- wavelength limit (q y — > 0), this evaluates to e scr ~ 
1 + (e 2 /27reoft) x Pgnr(£f) hi(2 /q y W), whereas in the 
metallic limit (q y W » 1), it is given by e scr w 1 + 
e 2 /7r 2 €oKhvF' The general form given in Equation 13 is 
used for the calculations that follow. 



III. SCATTERING MECHANISMS 

Low-field carrier transport in GNRs is affected by var- 
ious scattering mechanisms. Among them, scattering by 
acoustic and optical phonons, charged impurities, and 
edge-roughness scattering are expected to be most effec- 
tive. Other mechanisms such as remote optical phonons 
due to polar coupling with underlying substrates or di- 
electrics may be present, but are not considered in this 
work. Among the scattering mechanisms considered, 
phonon scattering sets the intrinsic limit on carrier mo- 
bilities. For GNRs with imperfect edges, line-edge rough- 
ness (LER) scattering can be rather strong for very thin 
ribbons. Unintentional charged impurities present either 
attached to the GNRs or spatially separated from them 
(for example embedded in the dielectric surrounding) also 
degrade the mobility. These scattering mechanisms are 
considered in this work; their effects are compared, and 
their relative importance is qualitatively and quantita- 
tively evaluated. 



A. Acoustic Phonons 

For acoustic phonon scattering, the zone-center acous- 
tic phonon interaction with carriers result in intra-valley 
scattering, which is quasi-elastic. Here we only consider 
the longitudinal mode since this mode induces higher de- 
formation potential than the out of plane and flexural 



modes [12]. The perturbation potential introduced by 
acoustic phonons is given by 
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where n = l/(exp (hco ac /kBT) — 1), n + 
v s q y is the acoustic phonon frequency, D, 
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n ,Uac = 

^ 16 eV is the 

deformation potential of acoustic phonons [11], p ~ 7.6 x 
10 -8 g/cm 2 is the 2D mass density of graphene, and v s ~ 
2 x 10 6 cm/s is the sound velocity in 2D graphene. Using 
the formalism outlined in section II leads to a scattering 
matrix element 



V ac (q y ) = D ac \q y \(l + e^') ] 
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where the Kronecker delta function S q ^k ensures 
the momentum conservation for the carrier + acoustic 
phonon system. The resulting intra-subband scattering 
rate is given by 
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where n p h = n + + n~ is used to take both absorption 
and emission of acoustic phonons into account. Here £ is 
the energy of carriers measured with respect to the Dirac 
point. The scattering rate may be re- written in the form 
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where backscattering restricts the value of \q y \ = 2\k y \. 
This form highlights the proportionality of the scatter- 
ing rate to the DOS of the GNR. The scattering rate is 
also inversely proportional to the width of the GNR, and 
becomes more severe for narrower ribbons, similar to the 
diameter dependence in carbon nanotubes [12]. 



B. Optical Phonons 

For optical phonon scattering, we consider the zone- 
boundary phonon of energy Hujlo ~ 160 meV, with an 
optical deformation potential D op ~ 1.4 x 10 9 eV/cm, 
assumed to be the same as for 2D graphene. Optical 
phonon scattering is inelastic; at low bias voltages and 
small electric fields, optical phonon emission by carriers 
is energetically forbidden, and absorption is damped by 
the high energy and low population of LO phonons. Since 
the emission of optical phonons is possible only if the 
kinetic energy of carriers exceeds Hulo, this mechanism 
is important only for highly energetic carriers. 

The perturbation potential introduced by LO phonons 
is given by 
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where n~ = l/(exp {h^Lo/^BT) — 1), n + = 1 + n~ , 
similar to those of acoustic phonons. The matrix element 
for optical phonon scattering is given by 



V op (q y ) = D op (l + e id ^'] 
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and the final scattering rate is given by 

° P Pgnr(£')Q + COS0 fefe ')> 



r op (£) ApWuo LO 



(19) 



(20) 



where the energy of the carrier after the scattering event 
changes to £' = £ ± fk^o • 



C. Line Edge Roughness Scattering 

To capture the effect of line-edge roughness of the 
GNR on charge transport, the width of ribbon is treated 
as a function of the longitudinal axis y. The width is 
given by W(y) = W + 5W(y), where SW(y) describes 
the roughness and W is the spatially averaged width 
(< W(y) >= W , and < SW(y) >= 0). This is shown 
schematically in Figure 1. The edge roughness SW(y) 
can effectively be described by two parameters and by 
an exponential spatial correlation function 
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where H is the amplitude, and A is the correlation length 
of the roughness. The LER leads to a spatially modu- 
lated bandgap, and the resulting fluctuations in the band- 
edge potential cause the scattering of carriers. The per- 
turbation potential for the nth subband is given by 
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(22) 



where £ n = hvp \k n \ is the conduction band energy of the 
nth subband relative to the Dirac point. 

Since the scattering potential is only dependent on the 
longitudinal axis, the scattering is intra-subband. The 
square of matrix element of the nth subband is given by 
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Since £ n ex 1/W, the scattering rate is proportional 
to 1/VK 4 , and if the LER scattering rate is dominant, 
the mobility should scale with the width of the GNR 
as Pler ^ W A . This behavior contrasts with inter- 
face/surface roughness scattering in traditional semicon- 
ductors with parabolic bandstructures. Since energy 
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FIG. 1: Schematic representation of line-edge roughness 
(LER) scattering in GNRs. The variation of the GNR width 
leads to spatial variation of the bandgap, and the fluctua- 
tions in the band edges cause carrier scattering and reduction 
of mobility. 
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FIG. 2: Line-edge roughness scattering rates in GNRs for 
(a) various roughness amplitudes, and (b) various roughness 
correlation lengths. 



eigenvalues in such traditional semiconductors are pro- 
portional to W~ 2 , confinement of carriers into a length 
scale of W leads to a roughness scattering limited mo- 
bility which scales as pm ~ W 6 [13]. Thus, the band- 
structure of GNRs make them inherently more robust to 
LER scattering than parabolic bandgap semiconductor 
nanostructures of comparable size. 

In addition to the the W~ 4 dependence of the LER 
scattering rate, the dependence on the roughness ampli- 
tude is H 2 . The factor A/(l + 4/c^A 2 ) is maximized for 
2k y ~ 1/A, indicating that LER scattering is the most 
severe for those carriers in the GNR that have Fermi 
wavelengths of the same order as the correlation length 
of the fluctuations. 

The LER scattering rates evaluated for realistic GNRs 
are plotted in Figure 2 (a) and (b). Figure 2 (a) shows the 
LER scattering rate as a function of the roughness am- 
plitude, ranging from ~1% to ~50% of the GNR width, 
for nominal GNR widths W = 3, 5, 8, & 12 nm. The 
correlation length for this plot is A ~ 3 nm. Figure 2 (b) 



illustrates the dependence of the scattering rate on the 
Fermi- level; as the Fermi level in the GNR is increased by 
a gate voltage, the Fermi- wavelength of carriers decrease, 
and therefore, LER scattering becomes more sensitive to 
roughness at shorter correlation lengths. Therefore, for 
2k y h « 1 and for 2k y h » 1, the LER scattering rate 
is low. 

We note here that for very rough GNR edges charac- 
terized by H ~ W ~ A, it is possible that carriers are 
localized in quantum-dot like confining potentials; trans- 
port in such rough GNRs will occur by hopping and tun- 
neling between localized states as opposed to diffusive 
band-transport (see [14]). The transport treated in this 
work is restricted to the diffusive band-transport regime. 



D. Impurity Scattering 

Due to an impurity charge located at (xq,0, z) with 
respect to the GNR which is on the z=0 plane (see Figure 
3), the unscreened scattering potential experienced by a 
mobile electron located at (x, y, 0) in the GNR is given 
by the Coulomb potential 
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where z is the distance of the impurity from the GNR 
plane and k is average relative dielectric constant of ma- 
terials on the two sides of the GNR [15]. The distance of 
the impurity from the origin is d = ^ z 2 + (x — xo) 2 . The 
screened impurity perturbation matrix element is then 
given by 



Vcoui(q y ) = (l + e iekk ') 
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Ko(\5k y \d)dx, 

(26) 
- k y f is the change in the carrier wavevec- 
tor along the GNR axis upon scattering, and KqQ is 
the zeroth order modified Bessel function, and e scr is the 
screening factor in Equation 13. 

In order to get the scattering rate due to all impuri- 
ties, we integrate over the distribution of impurities. For 
the most general distribution of impurities, the scattering 
rate is given by 
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where S is an effective scattering 'cross-section' with di- 
mensions of length, given by 
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(28) 
Here nsd(^) is the volume- density of the impurities which 
can vary with the distance from the GNR plane, t is the 
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FIG. 3: Schematic representation of the GNR and charged 
impurities that cause scattering of mobile carriers. 



thickness over which the impurities are distributed. This 
general formalism allows us to evaluate, at the same time, 
the effect of scattering by volume-distributed impurities 
as well as impurities located at the GNR surface. If the 
charged impurities are located at the GNR/dielectric in- 
terface or on the GNR surface, then the scattering rate 
is found unambiguously by taking z — > and nsd(^) — > 
^2dS(z), where ri2d is the 2D impurity density. For the 
evaluation of the scattering rates in the rest of the pa- 
per, we consider the volume (bulk) and surface impurities 
separately to highlight their relative importance. This is 
motivated by the recent finding that for 2D graphene, 
surface impurities are responsible for low-field mobility, 
and when they are removed, an order of magnitude im- 
provement in mobility is observed, even at room temper- 
ature [17, 18]. 



IV. CARRIER MOBILITY 

The rates for each type of scattering mechanism are 
first evaluated as a function of carrier energy to gauge 
their relative importance. They are plotted in Figure 4 
for a GNR of W = 5 nm width at T = 300 K. The 
relevant parameters used for the calculation are a edge 
roughness characterized by (H, A) = (0.5, 3) nm, a bulk 
impurity density nsd — 10 15 /cm 3 , and a surface impu- 
rity density ri2d = 10 10 /cm 2 . In addition to the individ- 
ual scattering rates, the total scattering rate obtained by 
Matheissens' rule as a sum of the individual scattering 
rates is also shown. For the specific case of the param- 
eters chosen, the LER scattering is seen to dominate at 
low carrier energies, followed by coustic phonon and sur- 
face impurity scattering. The effect of bulk impurities is 
found to be relatively weak compared to LER and sur- 
face impurities. It is possible to experimentally lower (or 
eliminate) scattering due to LER, bulk, and surface im- 
purities, as they are not intrinsic to the GNR, and such 
techniques are expected to be developed in due course. 

The intrinsic scattering mechanisms - due to acoustic 
and optical phonons, are also shown in the figure. Due 
to the dependence of the scattering rates on the ID DOS 
of the GNR, the general trend is a decrease of scatter- 
ing rates as the energy of the state increases from the 
band edge (£ c )> an d then the appearance of a step at 
the onset of the next subband. When the energy of a 
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FIG. 4: Scattering rates due to various scattering mechanisms 
for a 5 nm wide GNR as a function of the carrier energy. 
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FIG. 5: Carrier mobility in a 5 nm GNR as a function of 
temperature for three distinct Fermi levels. The relative im- 
portance of the various scattering rates can be ascertained 
from the plots. 



state is £ c + huiLO, optical phonon emission is allowed, 
and for such high energy states, optical phonon scattering 
dominates over all other scattering mechanisms. This is 
especially important at high-bias conditions, when high- 
energy states are occupied. Optical phonon scattering is 
responsible for the saturation of current flow through the 
GNR, leading to a sharp degradation of carrier mobility. 
The same has been found earlier for carbon nanotubes 
[16]. We restrict the rest of the discussion here to low- 
bias mobilities. 

Using the formalism for the calculation of carrier mo- 
bilities outlined in section II, the carrier mobility is calcu- 
lated for a 5 nm wide GNR as a function of temperature 
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FIG. 6: Carrier mobility in a 3 nm GNR as a function of 
temperature for three distinct Fermi levels. 



for three different locations of the Fermi level. The impu- 
rity densities used in this calculation are nsd = 10 15 /cm 3 
and ri2d = 10 10 /cm 2 . The results are shown in Figure 
5. When the Fermi level is located in the bandgap at 
£f = 0.75£ c from the Dirac point, the carrier mobil- 
ity is severely affected by surface impurity scattering at 
low temperatures as seen in Figure 5(a). Above ~200 
K, acoustic phonon scattering dominates the scattering 
rate. The dominance of surface impurity scattering for 
this case is due to the fact that when £p = 0.75£ c and 
the temperature is low, the net mobile electron density 
in the conduction band is low, and screening is weak. As 
the Fermi level is raised to the conduction band edge - 
£p = £ c (Figure 5(b)) and then above it - £p — 1.25£ c 
(Figure 5(c)), surface Coulomb impurities are strongly 
screened, and LER scattering dominates at low tem- 
peratures. At and around room temperature, acoustic 
phonon scattering still limits the electron mobility. The 
effect of LER scattering is felt only when the surface 
impurity density is low. The surface impurity scatter- 
ing rate scales inversely with the density of impurities; 
thus if the surface impurity density is much higher, say 
n2d ~ 10 12 /cm 2 , it becomes the dominant scattering 
mechanism, limiting the room temperature mobility to 
10 — 1000 cm 2 /V.s. Mobilities in this range has been 
recently reported in ultrathin GNRs [6], but due to the 
lack of experimental evidence of the GNR edge roughness 
of such samples, it is early to make a direct comparison. 
However, from the analysis, it can be concluded that such 
values of mobilities are not intrinsic, and result from ei- 
ther LER or surface impurity scattering. 

If the width of the GNR is lowered to W = 3 nm but 
with a roughness the same as the 5 nm wide GNR, LER 
scattering becomes stronger. Figure 6 depicts the relative 
importance of the various scattering mechanisms at dif- 
ferent temperature for three locations of the Fermi level. 
As can be seen, LER scattering starts dominating at and 
around room temperature as well, overtaking acoustic 
phonon scattering. If one considers the case when sur- 
face and bulk impurities are absent, the mobility will be 
determined by the relative importance of acoustic phonon 
and LER scattering. 
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FIG. 7: For very narrow ribbons, LER scattering is seen to be 
the mobility-limiting scattering mechanism, but wider GNRs 
- even for widths greater than ~5 nm are relatively insensitive 
to LER scattering. For such ribbons, the mobility at room 
temperature is limited by acoustic phonon scattering. 



Finally, if the surface impurity density is lowered to a 
level << 10 10 /cm 2 , the LER and acoustic phonon scat- 
tering mechanisms are of the most interest. Under such 
a condition, the mobility is plotted as a function of the 
GNR width for a constant edge roughness in Figure 7. 



The crossover of the LER and acoustic phonon scat- 
tering rates occurs at W ~ 4 nm for the chosen edge 
roughness parameters. For narrow ribbons, LER scat- 
tering dominates, but as the width of the GNR is in- 
creased, the LER scattering rate decreases as W A , mak- 
ing the GNRs of widths 5 nm and above relatively in- 
sensitive to edge roughness scattering. As W — > oo, the 
acoustic phonon scattering rate limit on the mobility ap- 
proaches that of 2D graphene; indeed, mobilities as high 
as 120,000 cm 2 /V.s have been recently observed in sus- 
pended 2D graphene sheets when the surface impurities 
were removed by current-induced annealing [17]. 



V. CONCLUSIONS & ACKNOWLEDGEMENTS 

In conclusion, it is shown from analytical modeling 
that low-field carrier transport in GNRs is affected by in- 
trinsic scattering by acoustic phonons, in addition to line 
edge roughness and impurity scattering. Impurities stick- 
ing to the GNR surface or at the interface of GNRs and 
underlying substrates are much more deleterious than 
those embedded in the underlying substrate. For very 
thin GNRs, the mobility is degraded by LER scatter- 
ing, which reduces the mobility as the fourth power of 
the GNR width. For many technologically relevant GNR 
widths, LER scattering is weaker than the intrinsic acous- 
tic phonon scattering. 
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